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1.1 Introduction 

After several decades of research in the development of a 
practical method for measuring acoustic intensity it was only 
recently that significant progress was made [1,2,3]. One, and 
probably the most important, reason behind this progress was the 
development of digital equipment for random data analysis capable 
of processing rapidly changing signals. In particular, the Past 
Fourier Transform Analyzer has made it possible quite easily to 
perform measurement of sound intensity. Another important rea- 
son, ha3 been the need to develop new methods of identifying 
noise sources and paths of propagation. The research has lately 
tf?ken three paths associated with three differed. methods in 
obtaining the intensity inside of and radiated from a structure 
namely by use of, the surface intensity, the acoustic intensity 
and structural intensity methods. 

All these methods are based on the estimate of the cross 
spectral density between two transducers and are therefore some- 
times referred to as cross spectral methods. For the case of 
acoustic intensity the cross spectral density is calculated 
between two closely spaced microphones, for surface intensity 
between a microphone and an accelerometer and for structural 
intensity between two closely spaced accelerometers. 

The surface intensity method [8-14] has been shown to be a 
useful tool in determining noise sources. One drawback is that 


this method requires measurements at many microphone / accoleronr 
eter positions which makes the measurements on a big structure 
somewhat time consuming* Another drawback is the problem associ- 
ated with the calibration of and phase shift between the micro- 
phone and accelerometer. However, a major advantage of the sur- 
face intensity method is that it gives the mean square surface 
velocity, radiated sound power and hence radiation efficiency of 
a vibrating surface all in one set of measurements. 

The acoustic intensity method [1-7] does not have the draw- 
backs mentioned. Unless safety considerations prohibit, the 
microphones can be swept over the surface giving simultaneously 
both a spatially and time averaged-signal. The calibration as 
well as phase shift are easy to account for in a measurement pro- 
gram. This method of measurement of acoustic intensity has also 
opened the way for several important applications in the field of 
acoustics, which previously required more expensive facilities. 
This was demonstrated in our work in 1980 when this new method 
was successfully used in the evaluation of the sound transmission 
loss of structures. [4,5] 

Very recently instrumentation has been developed by Bruel 
and Kjaer which is capable of evaluating and displaying in real 
time the acoustic intensity in one-third octave bands by using 
two closely spaced microphones. Another important factor, which 
has been carefully studied, is the calibration of and the ways of 
circumventing errors due to phase mismatch between the transduc- 
ers [6,7]. The estimate of the cross spectra obtained from an 


- 3 - 

PFT Analyzer may be subject to some statistical error e. Some of 
these errors are discussed in Reference 18, 

It seems that least amount of work has been done on the 
method of measurement of structural intensity, [15,16,17]. This 
method might be v'-.r/ useful in determining reflection and 
transmission coefficient 4 and coupling loss factors across struc- 
tural joints. 

Part one of this report presents both analytical and experi- 
mental results for the transmission loss of complex structures. 
On one side of the specimen the acoustic field is diffuse. On 
this side both the surface and acoustic intensity technique would 
give an erronous result when measuring the intensity since the 
pressure signal would then be composed of both incident and 
reflected soundwaves. On the other side of the specimen it was 
decided to use the acoustic intensity method. 

In this chapter a brief review of the literature on surface 
intensity, acoustic intensity and structural power flow is 
presented. 


1 , 2 Surf ace_ Intensity 

Some important progress [8-14] has been reported on both the 
theory and application of the surface intensity technique as a 
tool for determining noise sources. The equation for the inten- 
sity is [13] 
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1 " 2^ l (tl [°p« coa * + C p« aln *J dt <1-1) 

where Q pft is the imaginary part of the one-sided crose- spectral 
density between pressure and acceleration , C pa is the real part 
of the one-sided cross spectral density between pressure and 
acceleration and ♦ is the instrumentation phase shift* The main 
advantage with this technique is that it not only gives the radi- 
ated power but also the mean square surface velocity. Knowing 
both these quantities it is possible to calculate the radiation 
efficiency of the structure. It is important to note that a 
knowledge of the surface vibration is essential when attempting 
to reduce the emitted noise. One important drawback is that the 
method only works for the case of structure-radiated sound. 


1 • 3 Acoustic-Intensity 

It appears that the first successful attempt to use the 
cross spectrum between two microphones was made by Alfredson [1], 
He based his formulation assuming a simple harmonic signal. A 
more general expression was derived separately by Fahy [2] and 
Chung and Pope [3] . Probably the most useful expression for the 
acoustic intensity in connection with the use of a fast fourier 
analyzer was derived by Chung [6], The equation for the inten- 
sity is as follows: 


I - Im [ G 12 )/P w Ar 


( 1 . 2 ) 
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where |m ( ) stands for the "imaginary part of" and is the 
one-sided cross-power spectrum p is the density, u is the angu- 
lar frequency and Ar is the microphone spacing. Chung also 
showed that by measuring the intensity twice, where the micro- 
phones and connections to the FFT were switched in the second 
measurement, he could circumvent the error caused by instrument 
phase mie-match. Another method of intensity measurement has 
been developed where the phase and gain mis-match errors are 
eliminated. This is accomplished by a microphone calibration, 
and consequently the measurement time is reduced by about a fac- 
tor of two. 

The errors associated with the finite difference approxima- 
tion of the pressure gradient between the microphones have been 
carefully studied by Pavic. [15] 


1 , 4 Structural_Intensity 

Quite an extensive study of the measurement of intensity of 
waves in structures is given by Pavic [16]. He has formulated 
methods for the cases of one- and two-dimensional wave propaga- 
tion. The signal processing is done in the time domain with 
vibration transducers other than accelerometers. Similar work 
using accelerometers has been done by Verheij [17], in which the 
result is in the frequency domain. The equations for the power 
flow are expressed in terms of the spacing and cross power spec- 
tral density between two closely spaced accelerometers. The 
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power flow equations for bending waves on a beam are given by 


<P B>t 


2 (Bin* ) 
A 


2 - im jc ] 


df 


w 


( 1 . 3 ) 


and for longitudinal (L) and torsional (T) waves on a beam-like 
structure are, 

< P v - r§E " df (1.4 

L t A o 

<P v ~ zis / lm l G , lag (l.s 

where a^ and a 2 are the transverse accelerations, and and <* 2 

are the angular accelerations of the accelerometers, G(a^,a 2 ,f) 
is the one sided cross spectral density between the accelerome- 
ters, S is the cross sectional area, E is Young's modulus of 
elasticity, T is the torsional stiffness, B is the bending 

stiffness, m* is the mass per unit length and A is the spacing 
between transducers. The same reference also gives a table where 
equivalent signal processing operations in the time and frequency 
domain are presented. 
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2. THE APPLICATION OF THE ACOUSTIC INTENSITY TECHNIQUE TO THE 
SOUND TRANSMISSION LOSS OF PANEL STRUCTURES 

2.1 Introduction 

In this chapter a new approach for evaluating the transmis~ 
s ion lose of panels is presented. This approach is based on an 
estimate of the cross power spectrum between two closely spaced 
microphones [6]. The measured results will be compared with the 
transmission loss obtained by the conventional method), as well as 
the well known mass law. It will also be shown that by using 
this approach it is possible to identify paths of acoustic energy 
flow through a composite panel structure. 


2 , 2 Histor ical_Backqround 

The transmission loss is defined as ten times the common 
logarithm of the sound power incident on the structure divided by 
the transmitted sound power. In principle the receiving side 
should be anechoic [19]. It seems that Buckingham [20] presented 
the first theoretical formulation of the transmission loss of a 
partition separating two reverberant rooms. This method which is 
the current standard, requires two reverberant rooms. The 
transmission loss is given by the difference between the space- 
averaged sound pressure levels of the two rooms. However, there 
seems to be some scatter in the data obtained by this method even 
for laboratory measurements with a simple panel [21]. The same 
problem for the case of field measurements is discussed in 
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reference [22]. Besides these difficulties the method is also 
inapplicable whenever the receiving side is a free space or is 
very dead. Suggestions to overcome this problem are discussed in 
references [23,24]. Quite an extensive survey on available 
analytical models together with a comparison with measured data 
are presented in references [25,28]. 

The new method [4,5] to be discussed, has advantages over 
the conventional method since only one reverberant room is 
required and that corrections for the panel area and absorption 
on the receiving side are not needed. The new method also pro- 
vides a means of determining the sound power transmitted through 
different parts of a composite structure* In this report a sin- 
gle and a composite panel were chosen to verify the new method. 
A comparison was made between the measured transmission loss 
obtained by using the acoustic intensity technique and that 
obtained by the conventional method and the mass law. The sur- 
face intensity [13,14] technique would be an alternative way of 
measuring the transmitted sound power. 


2 . 3 Transmi88ion_Lo88_of_ a _Panel_and_it8_Mea.iurement 

The difference between the incident and transmitted inten- 
sity level is needed to rate the insulation of a panel [25]. 
There is no doubt that the transmission suite method is the most 
frequently used method for measurement of the transmission loss 
[27]. In this method the two reverberant rooms should be well 
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vibration-isolated to prevent mechanical tranemieeion between 
rooms. The volume of a reverberant room sets the lower limiting 
frequency for a diffuse field. The excitation in the transmis- 
sion room should be steady. By either using a rotating boom in 
each room or Bimply measuring simultaneously over an ensemble of 
microphones or by repeating the measurement for various micro- 
phone positions , space-time-averaged sound pressure levels in 

each room L_ and L. can be obtained. The noise reduction is 
P 1 p 2 

then defined as: 


NR - L - L * (2.1) 

P 1 p 2 

Provided that the flanking transmission can be neglected, the 
transmission loss is given by: [27] 


TL 


NR + 10 log 


A cT 
—R — 2 _ 


24 V 2 Ln(10) 


( 2 . 2 ) 


where A is the area of the panel, T 9 is the the reverberation 
time of the reception room, V 2 is the volume of the reception 
room and c is the speed of sound. The mass law is given by: 


TL - 20 log (rof ) -48 dB (2.3) 

where: m - mass/unit area (SI units), and f is the frequency, 

Hertz. 

q 

The reception room was simulated by a large wooden box of 2,8 m 
volume. Care was taken to seal the box. The sound pressure 
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level was averaged over 9 microphone locatione in the box. 

In the caee of a composite panel with several elements com** 
posed of different materials having different transmission coef~ 
ficients r., i-1,2,3... and corresponding areas A^, i«l # 2,3 f ... , 
the average transmission coefficient is given by [29] 


T 


j{ E Vi 


(2.4) 


where A p " E is the total panel area. The transmission loss 
of the composite panel is then: 



10 log 
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for a two element wall it can be shown [29] that the transmission 
loss is given by 
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where 
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2*’* Pi f C u» id^In^ tic—Intens ItV 

The average rate of flow of sound energy through and normal 
to a unit area in a sound field is by definition, the sound 
intensity in that direction [31] * The energy density in a dif- 
fuse field is given by [30] 


€ 


j2 

■ rms' 


p;« n/pc‘ 


(2,7) 


and the space-time averaged intensity is given by [30] 


<I> 




( 2 . 8 ) 


where p rmB is the rms sound pressure, p is the density of air and 
c is the speed of sound. 


2 . 5 Acoustic_IntenBity_by_Cross^.Spectral „ Method. 

Consider two omni-directional pressure microphones separated 
by a distance t , For a random, stationary, ergodic signal [31], 
the acoustic intensity I at a location in a sound field in a 
given direction r ia [5,6] 

t - E{pu t ) (2.9) 

where E[ } stands for the expected value and p and u are the 

JC 

acoustic pressure and particle velocity respectively. The pres- 
sure in equation (10) is the pressure at the midpoint between the 
microphones. It can be approximated by the mean value 
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Pi + Po 

p : -1— . ■ « (2.10) 

where P if i * li» 5 is the signal response from the microphones. 
By using a centered finite difference approximation [36] it can 
be shown that the particle velocity can be expressed as 


r * p Ar w 


( 2 . 11 ) 


where w is the angular frequency. This derivation is based on 
the criterion that fcAr << 1. Using equations (10), (11), and 
(12), it can be shown that the intensity is given by [6] 


Ij. * Im [o X2 ] / p A r tj (2.12) 

where Im( ) means the imaginary part of the argument, G^ 2 is the 
one sided cross power spectrum between the two pressure signals, 
fly using the transfer function between the two microphones, the 
correction for the phase shift can be included in the measurement 
method. The corrected cross power spectrum is [5] 



G 1 .2^jL2 


(2.13) 


where G^ is the measured cross spectrum between the two micro" 
phones, T^ 2 is the transfer function between microphone systems 1 
and 2 and |H,i l 0 the Gain of microphone system 2. For a 
description of the procedure in measuring the transfer function 
see the next section* 
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2.6 Expat iroental-Measurmnents 
2.6.1 Introduction 

The diffuse field in the reverberation room was measured by 
taking about 100 samples of the pressure signal from a rotating 
microphone boom during which time it made about two revolutions. 
The transmitted intensity was obtained by using the acoustic 
intensity method. The pair of microphones were swept with the 
center line of the pair held perpendicular to the panel surface. 
The distance between the midpoint of the microphones and the 
panel surface was about 20 mm. The speed with which the sweeping 
was done was approximately 4 em/see. The panel was clamped in an 
iron frame and fixed in the window of the reverberant room. The 
edges were sealed with clay to reduce the flanking transmission. 
A schematic diagram of the experimental set up is shown in Figure 
2.1. A diffuse field was used for the excitation of the single 
panel. This field was generated by an air jet noise source which 
gave most power ir. che higher frequencies and a loud speaker fed 
by a random noise generator which provided high levels for the 
lower frequencies. The spectrum was fairly flat and the overall 
level came to about 110-120 dB. During later measurements on the 
composite panel, the loudspeaker was replaced by a pneumatic 
driver which increased the levels at the lower frequencies and 
raised the overall sound pressure level (SPL) in the reverberant 
room to about 133 dB. 



Figure 2.1 Experimental set-up for measurement of trans- 
mission loss of a panel by acoustic intensity method. 
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2.6.2 Measur ement_of _Acoust ic_ I ntene ity 

The experimental set-up for measuring the transfer function 
between the two microphones is shown in Figure 2.2. The 1/2 inch 
microphones were calibrated before they were inserted at the end 
of the tube. A loud speaker driven by a random noise generator 
was fixed at the other end of the tube. The tube was made long 
enough to assure that plane waves existed at the location of the 
microphones. The signals from the microphones were fed through 
two measuring amplifiers before being processed by the FFT sys- 
tem. The transfer function was evaluated by using a key-board 
program. 

The experimental set up for measuring the acoustic intensity 
is shown in Figure 2.1. The array of two microphones were taped 
together. Their separation distance was about 14 mm as discussed 
in section 2.6.1. The acoustic intensity was then computed from 
the cross power spectrum of the two microphones as given in equa- 
tions (13) and (14). 


2.6.3 Measur ement_of_Free-Field_Intens ity 

The two microphone intensity technique was tested in the 

following way. The plane wave intensity was measured in an 

anechoic room by using the two-microphone technique. One hundred 

samples of the pressure signals were taken. This result was com- 

2 

pared with the plane wave intensity P rm8 /P c * evaluated from the 
auto-power spectrum from one of the microphones. 


The noise source was a loudspeaker fed by a random noise genera- 
tor. The measurements were made at a distance of 1,5 m from the 
loudspeaker. Figure 2.3 showB a comparison between the two 
methods in obtaining the intensity. This shows that the agree- 
ment is good between 400 and 2500 Hz. For frequencies below 400 
the discrepancy is about 0.5 dB and above 2500 Hz. it is about 
1.5 dB, 


2 . 7 Measurement of Diffuse Field Intensity 

A second test was performed to verify the two-microphone 

2 

method. The diffuse field intensity p„„/4pc, was measured 
inside the reverberation room while the window was open. The 
noise source consisted of a jet noise source run in combination 
with a loudspeaker. The auto-power spectrum of the pressure sig- 
nal measured by a microphone on the rotating boom was sampled 100 
times. The intensity at an imaginary surface in the open window 
in one wall of the reverberant room was then obtained by the 
acoustic intensity method. The signal was sampled at the same 
rate as in the previous measurement. 

Figure 2.4 shows the plots of diffuse field intensity in 
narrow bands obtained by these two methods. Curve number 2 has 
been lowered by 10 dB to show clearly the very good similarity 
between the two curves. The measurements were made by 1/4 inch 
microphones. The corresponding plot in one-third octave band 
averaged values is shown in Figure 2,5. The discrepancy above 
4000 Hz. is most likely with the two-microphone results and is 
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due to the fact that kAr << 1 is no longer satisfied. Because of 
this , the pressure gradient is not accurately determined by the 
centered finite difference approximation (Equation 12). 


2.7.1 Measurement_of_.Tr anomies ion_LQBS 


Two different panel structures were selected for the 
transmission loss measurements. The first structure was a 
clamped steel panel of dimensions 121 cm x 121 cm x 1.6 mm. The 
area of the panel on which the sweeping was done, ie. the area 
facing the outside of the reverberation room excluding the part 
covered by the wooden frame was for the steel panel and all the 
other panels was 113 cm x 113 cm. The second structure was an 
aluminum panel with a thickness of 3.2 mm. The panels were 
clamped in a window of the reverberation room. The room volume 
was 212 cu.m. The sound pressure signal from the microphone on a 
rotating boom inside the reverberation room was used to determine 
the incident intensity. The transmitted intensity was obtained 
by sweeping the microphones on the outside of the panel. The 
result in one third octave bands is shown in Figures 2.6 and 2.7. 

The transmission loss of a composite panel with the same 
inside dimensions and 3.2 mm thick was also measured. The compo- 
site panel consisted of an aluminum panel with a plexiglass win- 
dow which covered about one ninth of the total panel area. The 
plexiglass window was concentric with this panel. The plexiglass 
thickness was about 1.6 mm. The result from the measurement of 







X 
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the transmission loes of the two different mater lain of the com- 
posite panel is shown in Figure 2.8. Also in the same figure are 
shown the predicted values as obtained from Equation (4). It is 
clear that this method is, able to identify the power transmitted 
through different areas of a complex structure. 


2.7.2 Results_and_Conclus ions 

The results obtained by using the acoustic intensity tech- 
nique agrees well with the transmission loss measurements made by 
the conventional transmission suite method as seen in Figures 2.6 
through 2.8. 

There are several advantages with this new method Firstly, 
it only requires one reverberant room. This means that the cost 
in building facilities for measuring transmission loss can be cut 
substantially. This also means that for already existing facili- 
ties with only one reverberation room there is now a method 
available for measuring the transmission loss. Secondly, it 
seems quite possible, with this method, to identify paths of 
power flow through structures. The conventional method does not 
have this possibility. When the measured values of transmission 
loss (by this new method) of the plexiglass and the aluminum 
parts of the compostie panel were substituted into Equations (4) 
and (5) the resulting total transmission loss agreed fairly well 
with the total transmission loss obtained by measurements made 
over the whole area of the composite panel. 




One of the problems associated with this method is the fre- 
quency limit set by the microphone spacing. This could be over- 
come by reducing the spacing to obtain intensity values in higher 
frequency ranges or increase the spacing of microphones to obtain 
data in lower frequency range. 


3. EXPERIMENTAL AND THEORETICAL STUDIES OF THE SOUND TRANSMIS- 
SION THROUGH SINGLE , COMPOSITE AND STIFFENED PANELS 

3.1 Introduction 

After the two microphone intensity technique had proved to 
be successful in measuring the transmission loss of a homogeneous 
single panel, it was decided that the same method should be 
applied to more complicated structures such as a composite panel 
consisting of a metal panel with a plexiglass window and a panel 
with stiffeners. The latter structures are more similar to those 
of a typical small aircraft fuselage than the homogeneous single 
panels studied earlier. It was also decided that the measured 
values should be compared with those obtained theoretically. The 
theoretical model was based on the principles of Statistical 
Energy Analysis (SEA). 

3 . 2 Exper imental_Work 

The instrumentation and the procedure used were quite simi- 
lar to those outlined in Section 2. However there were some 


27 


diff erences. 

A powerful driver shown in Fig. 3.1 was used instead of a 
loudspeaker. The sound sources shown in Figure 3.1 and 3.2 were 
run simultaneously. A typical sound source spectrum and level 
measured in the reverberation room are shown in Figure 3.7. The 
overall sound pressure level was about 133 dB. For the composite 
panel, to obtain the acoustic intensity transmitted microphone 
sweeping was performed with a distance of about 20 mm from the 
midpoint between the microphones to the panel (see Fig. 3.3 and 
3.4). The clamped stiffened panel is shown in figures 3.5 and 
3.6. The sweeping was done (see Figure 3.6) over a grid composed 
of very thin threads. The grid was used to assure uniformity in 
sweeping from the structure. However, in this procedure the 
sweeping was performed at a distance equal to the beam width 
plus aoout 10 mm, (between the microphone mid-point and the 
beams,) as shown in the figure below. 



^8 

Microphone, 

mid-point 











Figure 3.7 Typical sound pressure level spectrum inside the reverberation 
room due to the pneumatic driver (20 psia/6.0V) and the jet noise (80 psia) 
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The reason for this was that if the sweeping was done closer to 
the structure, the effects of standing waves in between the beams 
and scatter from the beams may affect the measurements. 


3 . 3 SEAJModellinq_of_the_Panels 

3.3.1 Per ivat ion_of_the_Tr anemias ion_LosB 

Consider the composite panel clamped in the window of the 
reverberation room. In this case the composite panel consists of 
two Bub panels namely the plexiglass window and the aluminum 
panel. The diffuse field in the reverberation room is meant here 
to model the exterior noise field of a fuselage, and the panel a 
possible substructure. Before applying the statistical energy 
analysis to the following problem, it should be said that this 
theory is based on the fact that in each bandwidth of interest in 
each subsystem there must be enough resonant modes present to 
make the concept statistically meaningful. Following references 
27 and 32 it seems logical to divide the problem into three sub - 
structures: 

1. The reverberation room, 

2. The plexiglass window and 

3. The aluminum panel, 

as shown in Figure 3.8. Conservation of energy in each subsystem 
requires that 


Fiaure 3 ." Block diagram representing 
[hree coupled systems of the composite 


power flow between 
panel . 
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i, 


n i,in " n i,diss * ^ n ij 


( 3 . 1 ) 


where the prime on the eum means that j«i is excluded; further 
the power flows in equation (1) are as follows; 


n i,in " ^ n P u ^ P ow «* the ith subsystem , and 


n i,diss " u v i E i, tot 


(3.2) 


is the dissipated power of the ith subsystem, where Vi is the 
internal loss factor and is the total energy, and finally 


"ij - ° "u E i - “ "ji E j 


(3.3) 


is the power flow between subsystem i and j, is the coupling 

loss factor. 


The number of resonant modes N^N.. in a bandwidth of Af in 
subsystems i and j and the coupling loss factors and must 
satisfy the consistency relationship: 


N i ”ij ■ N j ”11 


(3.4) 


Using equations (3 . 1) , (3 .2) , (3 . 3) , and (3.4), the power balance 
equations of the three system model become: 
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where f noni:e8 stands for the non**resonant power flow and and 
n 3,rad *«P*«»®nts hers the power radiated to an infinite half 
space. The non-resonant power flow is given by 


n nonres “ r nonres n inc 


(3.6) 


where f nonEe8 io the non-resonant transmission coefficient and 
n inc is the incident sound sound power on the side of the panel 
facing the reverberation room. Now it holds that 


n ■ r n * t a - rt tys, . t 

nonres nonres inc nonres p 4pC 


(3.7) 
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where <p 2 > t * s the ®P ace time averaged sound pressure, A p is the 
panel area, p is the density of the ambient medium, C is the 
speed of sound and V is the volume of the room. With this in 
mind, equations (3.5) can be written as: 
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or in matrix form as: 


(3.9) 
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The non-resonant transmission coefficient r , of each aub 

nonres/ 

area ia found by U8ing the mass law [33] 


10 log ~ 


x, nonres 


- 10 log {^] 


(3.10) 
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then an average tr anemias ion coefficient ia given by 
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(3.11) 


where A, are subareas and A ia the total panel area. The 
1 P 

transmission loss ia defined as 


TL - 10 log 
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inc 


n rad + n nonrea 


(3.12) 


for the composite panel this becomes 


l n 2,rad + n 3,rad + n nonresJ 


The radiated power is [27) 
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where for ith subpanel area, M_ is the mass of the panel, , 

P< i t r ao 


is the radiation resistance to half space. The derivation of 

k 2 n 


R r rad given in [34] . The corrected equation ia given in 
[33] . Substituting equation (12) for i-1,2, and Equation (7) into 
Equation (11) gives: 
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A_C 

3v~ E x 


.277 


.277 


‘2, rad 


M. 


e 2 + e 3 


P3 


A 


r nonres E 1 


(3.15) 


rearranging this equation gives the final expression for the 
transmission loss as 
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from Equation (8). This Equation shows that foe the case o f a 
single panel Equation (14) reduces to 


TL - 10 log 


nonres 


E 

_4V _iad 


V s 


M. 


(3.17) 


Two computer programs have been written to solve for the 
transmission loss given by Equations (14) and (15). A reduction 
scheme [37] was written to solve the set of linear equations. 


In order to solve for the energies in Equation (8), the loss 
factors would have to be known either through measurement or 
based on empirical data. The following section discusses these 
parameters and how they were obtained, for the particular set-up. 


3.3.2 Loss Factors 


The loss factors are the internal loss factor, coupling loss 
factor and radiation loss factor. The internal loss factor of 
the reverberation room is related to the reverberation time by: 

V - (3.18) 

where T is the reverberation time, and f is the center frequency. 

An automated program developed by flruel and Kjaer [36] was 
used, to measure the decay of the sound pressure level in the 
room and to calculate the reverberation time. The triggering of 
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the noieo source was controlled by the program through the inter- 
face between the source and a HP 9825A desk top calculator, The 
levels were calculated from a real-time analyzer BK 2131 and the 
data was then transferred to the calculator, (See figure 3,9.) 
The rotating boom had a speed of about 1 revolution per 16 
seconds. Four recordings, each based on 65 spectra spaced 88 ms 
apart were taken to cover one revolution. 

The internal loss factors of the single panel were obtained 
by using a program developed by Simon Wang (details are given in 
Part II of this report). The panel, mounted in the window to 
assure the same boundary conditions as those present in the meas- 
urement of the transmission loss was excited by a shaker. Five 
accelerometer positions were used to record the decay of the sig- 
nal. For each accelerometer position, five decay curves were 
taken in each one-third octave band. The set-up and the results 
for the aluminum panel are shown in figures 3,10 and 3.11. As 
seen from figure 3,11, the internal loss factor is much higher 
at the critical coincidence frequency (3987 Hz). The reason is 
that the above program gives the sum of the radiation loss factor 
and internal loss factor. Other investigators [27,37] have shown 
that it is probably only near the critical coincidence frequency 
that the radiation loss factor may exceed the internal loss fac- 
tor . 

It is planned to measure the internal loss factors for the 
other panels later this spring. The radiation loss factor is 
obtained from theory [27]. 








Fig. 3.10 Experimental setup for the measurement of the 
the loss factor of the panel 
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The coupling loss factor between the plexiglaae and the 
aluminum panels has been derived. However the derivation is 
similar to the given in reference [38]. The derivation is 
lengthy and will not be presented here. 

3.3.3 The Single Panel 

It can be seen from Figure 3.12 that the SEA-model predicts 
the transmission loss fairly well for frequencies up to the coin" 
cidence dip. The disagreement above the critical coincidence 
frequency might be partially explained by the fact that the 
microphone spacing (14 mm) no longer satisfies the condition that 
kAr < < 1 . However* other experiments which show the effect of the 

spacing on intensity measurements show in contrast that this 
error should be small. Another more likely reason is that the 
panel is radiating less efficiently than that predicted by the 
theory. An average value of radiation loss in each one third 
octave band was taken while evaluating the transmission loss 
shown in the figure. 

3.3.4 The Composite Panel 

The narrow band values of the sound pressure levels and 
intensity levels are shown in figures 3.13 to 3.19. It can be 
seen in Figure 3.13 that the intensty level is fairly smooth up 
to about 5,800 Hz. The ripples above that frequency are due to } 

the effect of crossmodes in the tube where the transfer function 
was measured. So, this effect is carried through the calculation 
to the final value of the imaginary part of the cross spectrum. 
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bands and expanded in the frequency range 0 to 2000 Hz 





expanded in the frequency range 0 to 2000 Rz . 



Figure 3.18 Sound pressure level, L p , and intensity level 
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The intensity level has a minimum around 2000 Hz and 
increases up to about 5000 Hz. For the same frequencies, the 
intensity level, shown in Figure 3.14, inside the reverberation 
room drops up to about 4000 Hz, after which it recovers. This 
demonstrates experimentally the effect of the transmission loss 
dip at the critical coincidence frequency which for the aluminum 
panel is 3987 Hz. In figure 3.15 the intensity level and the 
sound pressure level are superimposed. 

The intensity from the plexiglass area of the composite 
panel iB shown in figures 3.16 and 3.17. For frequencies around 
3500 Hz to 5500 Hz, it seems that less intensity, than might be 
expected, ' <» emitted from the plexiglass area. A possible reason 
for this is that the aluminum panel area absorbs energy from the 
plexiglass area. 

The intensity obtained for the whole area is presented in 
plot 3 3.18 to 3.19. The general trend as shown in Figure 3.18 
seems to follow that of the aluminum sub panel area (see Figure 

3.14) . 

The comparison between predicted and measured values is 
shown in Figure 3.20. There is fairly good agreement for fre- 
quencies between 163 Hz and about 4000 Hz. The reason for the 
discrepancy above the coincidence dip is possibly due to uncer- 
tainties in the estimates of the internal loss factor and radia- 
tion efficiency. 
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3.3.5 The Stiffened Panel 

It was decided that the first attempt to formulate the 
powerbalance equations should only consider the subareas shown in 
Figure 3.21, The power balance equation is shown in Figure 3.22. 
However, for the measurements the panel was divided into five 
subareas; 

a) the small plexiglass window (1), 

b) the area in between the windows (2) , (3) , (4) f £ (5), 

c) the big plexiglass window (6), 

d) the left stiffened bay (7) ,(8), and half of (9), 

e) the right stiffened bay (10), (11), and half of (9), 

where the bracketed numbers indicate corresponding areas in Fig- 
ure 3.21, The narrow band values for the soundpressure and 
intensity levels between 0 and 2000 Hz for different subareas are 
shown in figures 3.23 to 3.26. Figures 3.27 to 3.31 show the 
consistency between the three measurements. If these graphs are 
overlayed, it can be seen that they are all quite close to each 
other above 2000 riz. This suggests that there might be some con- 
tamination from neighboring areas. This might be due to measure- 
ments made with too big a distance between the microphones and 
the panel. Below 2000 Hz it seems that the stiffened areas give a 
higher transmission loss. 
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Figure 3.28 Comparison between three measurements of the transmission loss 
of the area in between the windows in the stiffened panel. 
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The comparison between theory and measurement is shown in 
figure 3.32. The internal loss factor of the plexiglass, was 
assumed as 0.1. The critical coincidence frequency of the alumi- 
num panel is about 9900 Hz. It is seen that the agreement 
between measured and predicted values is good between 100 Hz and 
5000 Hz. Outside this range the intensity values obtained with 
the one-half inch microphones taped together, are not valid. 


3 . 4 Comparison Between the Single, Composite and Stiffened Panel 

The transmission losses of the aluminum sub-area on the com- 
posite and stiffened panels are compared with those of the alumi- 
num panel in Figure 3.33, The aluminum part of the composite 
panel and the whole panel is 3.18 mm thick. This gives a criti- 
cal coincidence frequency of 3987 Hz. The corresponding values 
for the aluminum part of the stiffened panel are 1.28 mm and 9900 
Hz. Since, in the region below the critical coincidence fre- 
quency, the transmitted intensity is approximately halved if the 
thickness is halved for the same material (see Eq. 9), it follows 
that the graph corresponding to the stiffened panel should be 
lower, in fact it should be about 8dB lower which is the case 
for most of the one-third octave bands. As can be seen in figure 
33, the higher critical coincidence frequency is the reason why 
the coincidence dip has been shifted to higher frequencies for 
the stiffened panel. If these reasons are taken into account it 
seems that the three graphs agree fairly well. 
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A similar comparison, but for the plexiglass panel is shown 
in Figure 3.34. It seems that the agreement is fairly good for 
frequencies above 1000 H2. Finally, a comparison of the 
transmission loss of the three panels is shown in Figure 3.35. 
If the graph of the stiffened panel is raised by about 8 dB for 
the aforementioned reasons it agrees fairly well with the graph 
of the aluminum panel. It seems that the effect of the plexi- 
glass window in lowering the transmission loss is more pronounced 
for the composite panel. 


4. the prediction of the sound pressure level in a rectangular 

CAVITY CONNECTED TO A RADIATING PANEL 
4.1 Introduction 

The estimate of the sound pressure level inside cubic or 
rectangular cavities is not a new and unsolved problem. It has 
been discussed in several texts and papers [42,43,44]. Usually, 
expressions for the sound pressure level at an arbitrary position 
due to an arbitrary source distribution are determined directly, 
by using an expansion in eigen functions [43] and employing a 
green's function technique. The solution for cases, where the 
walls are acoustically hard, are relatively simple to obtain. 
However, for a finite wall impedance the eigen functions become 
complex and the approximate solution is difficult to obtain. 

Another approach is the image source theory [41]. In this 
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theory the reflections from the walls (which could have a finite 
and the room equation do not give a relation between the sound 
wall impedance) are modeled as image sources. Adding up all the 
contributions from all the sources then gives the sound pressure 
spectrum at an arbitrary position.. The equation is similar to 
the known room equation and is easy to program on the computer . 
The image source theory has been used to estimate the sound pres- 
sure level along the center line of a fuselage cavity in refer- 
ence [40]. No comparisons between theory and measurement were 
shown in it. The image source theory has been used in this pro- 
ject. The equation given in reference [41] is basically of the 
same form as that used in [40]. The equation has been extended 
here for a three dimensional case. This theory and a discussion 
between predicted and measured values are presented here. 

At this point it should be said that both the image theory 
and the room equation do not give a relation between the sound 
pressure level at a point and the spatial distribution of the 
impedance on the boundary of the enclosure. In both these equa- 
tions, the room constant only accounts for an average absorption 
coefficient which is based on the reverberation time. 

4 • 2 The Image Source Theory 

consider first a source located between four parallel walls 
(forming a corridor) as shown in Figure 4.1. S and OB stand for 
the source and observer respectively. The solid lines are direct 


or reflected sound paths, (from the source to the observer). The 
dashed lines indicate the location of the image sources. If j is 
an index in the Z direction, it is seen in Figure 1 that the Z 
coordinates for sources on the positive side of the Z axis are: 

Z* - Z+_ L + 2( - q j ) Z q + Zq.(H - zj j-1 , 2 , 3 . . . (4 . 1 \ 

where : 

q. i - j[l -(-l) 1 ] (4.2) 

Similarly the Z coordinates for the image sources on the negative 
side are given by: 


z. * * Zq.Z 0 - ^ j ] ^ ^ * ’ * (4.3) 

where H is the height of the corridor and Z Q is the real Z coor- 
dinate. The corresponding equations in the Y-direction are 

Y i “ Y i-1 + 2 ( 1_q i] Y o + Zq i (w-Y 0 ),i«l,2,3... (4.4) 

and 

Y i - *1-1 - 2(3 i y o - 2 Hi) Ho] 1 - 1 ' 2 ' 3 -" (4 - 5) 

where W is the width of the corridor and Y Q is the Y-coordinate 
of the real source and 

q i - | |i-(-l) i ] (4.6) 

If now reflections are allowed in the X-Direction, the corridor 
becomes an enclosure and two more equations need to be added: 
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*k ' X k-1 + 2 ( 1 -’k] X o * 2 ’k ( H_Z o] k " 1 ' 2 ' 3, • • <«•?> 

and 

X~ " X k-1 ~ 2q k X o " 2 (l" q JcJ( L " x 0 ]k-l,2,3.., (4.8) 

where L is the length of the enclosure/ X Q is the X-coordinate of 
the real source and 

q k - (4.9) 

The total sound pressure level Lp at the observer position OB 
(X,Y,Z) is given by the room equation: 



u 


+ 10 log 




(4.10) 


where is the sound power level caused by the source, Q is the 
directivity factor, R T is the room constant and r is the distance 
between the observer and the real source. In the image theory 
the room constant is given by 


n ijk 
if 


T ijk IMtt r^ k 
where 


n. .. 


pos . sources 


ijk n4ir X, 


ijk 


(4.11) 


neq. sources 



(4.12) 


n is the total acoustic power radiated by a real source, 
is the power emitted by the (i,j,k)th image source, 
given by 


and n ijk 
"ijk iB 
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(4.13) 

n ijk " n F(« ? ,k)*P(fl< 1 ,-k)*P(o u ,j)*F(a d# -j)*P(a ff -i)*P(a ji ,i) 
with 

F(a,h)-(l-a) 9(h) (4.14) 

and 


9(h) 


/ 



0 K » 0 , 

| h/2 | ; h even , 
integer | l/2h|+ 
^integer | l/2h| 


1 for h; odd; positive, 
for h; odd; negative, 


(4.15) 


where the subscripts r, 1, u, d, f and a stand for right, left, 
up, down, front and aft, a is the absorption coefficient given by 


. 60V 
• CST 


(4.16) 


where V is the enclosure volume, S is the interior surface area, 
T is the reverberation time of the enclosure and C is the speed 
of sound. 


The image theory is based on the following assumptions: 

1. each sound source is a simple point source with a constant 
acoustic power output, 

2. the acoustic response of each wall can be described in terms 
of an average sabine absorption coefficient , 

3. the plane walls can be replaced by a set of gradually fading 
image s, and 
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4, the contribution of all sources (real and image) are added 
incoherently , 

4 . 3 Interior Noise Level as a Function of Space and Frequency 

The image theory* discussed in the previous section* is now 
applied to the fuselage in the following suggested way: consider 

first the experimental set-up shown in figure 4.2, The idea here 
is that the structure is divided into substructures or more 
correctly into sub areas* It was decided that the division into 
subareas should be based on the structural differences within the 
panel. The plexiglass elements, the stiffened part of the panel 
in between the windows, as well as the area below the windows as 
seen in Figue 4.2 formed the different subareas. Since all the 
sources are considered to be incoherent* the room equation is now 
modified to account for a summation from all sources, 


' P ( r ijk' £ ) 


t 

n 


H.( r ijk' £ J + 10 109 


477 r 


ijk 


4 

vn 


(4. 17) 


where n is the number of subareas and ^(r^^f), is the radiated 
intensity from the n th element times the area. 


Figure 4.1 The location of observer (OB) and real source 
(S) position and image source locations, for the image 
source model. 
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4.4 Experimental Procedures 

The aide of the panel, in Figure 4.2, which ie facing the 
reverberation room was excited by a reverberant sound field with 
an overall level of about 130 dB. This level was generated by an 
electro-pneumatic driver and a jet noise source. For all meas- 
urements; the auto power spectrum from the rotating boom in the 
reverberation room and the signals from the microphone in the 
cavity as well as the cross power spectrum from the two micro- 
phones on the cavity side, were sampled 100 times. These spectra 
were then stored on the disc and later processed in a fortran 
routine based on equations (4.10) and (4.11). 

The pressure supplied to the electro-pneumatic driver could, 
for the given sound pressure level, only be maintained at a con- 
stant level with fluctuations of 10% for about five minutes. 
This had the result that the noise source had to be turned off 
and on between all measurements of the intensity level from the 
subareas, the intensity level in the reverberation room and the 
sound pressure level in the cavity. It is difficult to estimate 
the error inherited in this procedure. It is believed to be 
within IdB, for most frequencies. The sweeping of the two 
microphones was carried out as described in Chapter 3. 

4 . 5 Results and Discussion 


The results obtained by using the image theory as well as 
those obtained by the room equation were calculated and compared. 
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The room constant used was; 

4 ' -E S n ^ (l-« ) <♦•«) 

n 

where S n - eubarea and <* n ia the aborption coefficient from the n 
th sub area (diffuse field theory). In the following pages a 
preliminary output from the fortran subroutine for the case of 
the stiffened panel is given. In the output the following abbre' 
viations are used: 

SPLT « sound pressure level at the observer 
position (measured) 

IMTH ■ sound pressure level at the observer 

position predicted by the image source 
theory. 

DFTH - sound pressure level at the observer position 
predicted by the diffuse field theory 

INLV » intensity level from the sub area (source no. X) 

The last page of the out put shows the contribution from all 
sources. The number of sources (5 real + all images) was about 
30,000. The image theory and diffuse field theory for the 
evaluation of sound pressure level inside the cavity seem to 
agree fairly well with the measured values between 400 Hz and 
2500 Hz. 
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OUTPUT 

Table 1 

FROM THE IMAGE SOURCE PROGRAM 


source-co 
,13 .94 

room-dim 

1.78 1.27 1.22 .80 

obs-c 

.53 

sum-1 im 
15 15 15 

absorption-coe£f . 

.10 .10 .10 .10 .10 

.10 

2.20 

area 

2.26 2.26 2.20 1.42 1.42 



Contribution From Source No. 1 


FFe^.'-'irgTH 

— .... ..... .. - 

—tjFCT — 

1 

TWTH r 

........ t 

jjpjrfj 1~ 

_ 1 

INLV 

11 

163 

1 

85.46 | 

! 

75.00 ! 

1 

75.66 J 

81.25 

200 

100.00 

88.35 

89.01 ! 

94.60 

250 

95.87 

90.75 

91.41 | 

97.00 

315 

95.11 

87.32 

87.98 j 

93.57 

400 

95.28 

85.91 

86.58 

92.16 

500 

92.83 

84.71 

85.37 

90.96 

630 

93.34 

78.37 

79.03 

84.62 

800 

85.14 

71.48 

72.15 

77.73 

1000 

80.39 

70.63 

71.29 

76.88 

1250 

77.02 

65.16 

65.82 

71.41 

1600 

71.62 

60.27 

60.93 

66.52 

2000 | 

70.01 

57.52 

58.18 

63,77 

2500 

69.52 

56,54 

57.20 

62.79 

3150 

67.24 

52.61 

53.27 

58.86 

4000 

64.88 

46.92 

47.58 

53.17 

5000 

65.85 

43.07 

43.73 

49.32 

6300 

66.12 

53.93 

54.59 

60.18 

8000 

1 64.20 1 

48.98 1 

49.64 

55.23 
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Contribution From Source No. 2 




IMTri T 

’’WTO T 

l 

INLV 

163 

85.46 

74.80 

— "■■■■■■” ■ ■ t 

75.47 ! 

77.52 

200 

100.00 

93.89 

94.56 

96.61 

250 

95.87 

94.38 

95.05 

97.10 

3.15 

95.11 

88.89 

89.55 

91.61 

400 

95.28 

87.30 

87.98 

90. 02 

500 

92.83 

84.17 

84.84 

86.89 

630 

93.34 

80.11 

80.78 

82.84 

800 

! 85.14 

76.95 

77.62 

79.68 

1000 

80.39 

71.78 

72.45 

74,51 

1250 

77.02 

68.48 ! 

69.15 

71.20 

1600 

71.62 

64.23 

64.89 

66.95 

2000 

70.01 

61.15 

61.82 

63.87 

2500 

69.52 

61.08 

61.75 

63.80 

3150 

67 . 24 

57.22 

57.89 

59.94 

4000 

1 64.88 

51.68 

52.34 

54.40 

5000 

1 65.85 

49.91 

50.58 

52.63 

6300 

! 66.12 

57.49 

58.16 

60.21 

8000 

1 64.20 1 

53.17 1 

53.84 | 

55.89 


Contribution 

From Source No. 3 



Freq. Hz I 

i 

— splt — r 
. 1 


DFTH T 

1 . 

INLV 

"1 

163 | 

85,46 

. — 1 — 

78.06 ] 

■ u ■ . >| 

78.76 

81.77 

200 

100.00 

90.32 

91.03 

94.04 

250 

95.87 

94.41 

95.11 j 

98.13 

315 

95.11 

89.22 

89.92 

92.94 

400 

95.28 

90.70 

91.40 

94.42 

500 

92.83 

86.62 

87.33 

90.34 

630 

93.34 

83.75 

84.46 

87.47 

800 

85.14 

76.65 

77.35 

80.37 

1000 

80.39 

71.51 

72.21 

75.22 

1250 

77.02 

67.78 

68.48 

71.50 

1600 ; 

71.62 

62.88 

63.58 

66.59 

2000 

70.01 

60.41 

61.11 

64.12 

2500 ! 

69.52 

59.30 

60.00 

63.01 

3150 ! 

67.24 

57.22 

57.93 

60.94 

4000 

64.88 

51.53 

52.24 

55.25 

5000 

65.85 

51.77 

52.47 

55.48 

6300 

66.12 

56.33 

57.04 

60.05 

8000 1 

64.20 1 

49.80 | 

50.50 1 

53.51 






Contribution 

From Source No. 4 


Preq. Hz. 

SPLT 

T 

IMTH 1 

L_ 

DPTH 1 
U 

INLV 

163 

85.46 

! 

i 

1 

81.04 | 

1 

81.74 | 

80.59 

200 

100.00 

i 

95.49 

96.19 ! 

95.05 

250 

96.87 

1 

1 

i 

97.15 

97.85 

96.70 

315 

95.11 

92.47 

93.17 

92.03 

400 

95.28 

1 

i 

92.24 

92.94 

91.80 

500 

92.83 

1 

1 

87.81 

88.51 

87.36 

630 

93.34 


85.52 

84.38 


85.14 

77.46 

78.17 

77.02 

1000 

80.39 

I 

| 

75.18 

75.88 

74.74 

1250 

77.02 

1 

* 

69.34 

70.04 

68.90 

1600 

71.62 

! 

i 

i 

64.44 

65.14 

64.00 


70.01 


64.33 

63.18 

2500 

69.52 

I 

i 

63.61 

64.31 

63.17 

3150 

67.24 

i 

1 

57.58 

58.28 

57.14 


64.88 

1 

i 

52.08 

52.79 

51.64 


65.85 

1 

i 

51.90 

52.60 

51.46 


66.12 

! 

59.58 

60.29 

59.14 

8000 

| 64 . 20 

| 

55.45 1 

56.15 1 

55.01 


Contribution Prom Source No* 5 


Preq. Hz. T 

SPLT | 
i 

IMTH 1 

- 1 

DPTH T 
\- 

INLV 

— 

163 ! 

85.46 

84.96 

85.69 ! 

84.55 

200 ! 

100.00 

95.95 

96.68 

95.54 

250 

95.87 

96.70 

97.43 

96.29 

315 ! 

95.11 

93.86 

94.59 

93.45 

400 

95.28 

87.48 

88.21 ! 

87.07 

500 

92.83 

86.58 

87.31 

86.17 

630 

93.34 

83.53 

84.26 

83.12 

800 

85.14 

76.91 

77.64 

76.50 

1000 

80.39 

72.78 

73.51 

72.37 

1250 

77.02 ! 

69.31 

70.04 

68.91 

1600 

71.62 ; 

64.36 

65.09 

63.95 

2000 

70.01 

63.59 

64.32 

63.18 

2500 

69.52 

64.11 

64.84 

63.70 

3150 

67.24 

59.24 

59.97 

58.83 

4000 

64.88 

53.79 

54.52 

53.38 

5000 ! 

65.85 

54.18 

54.91 i 

53.77 

6300 | 

66.12 

59.74 

60.47 

59.33 

8000 1 

64.20 | 

55.74 | 

56.47 1 

55.33 
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Contribution Of All Sources 


Freq. Hz. 

1 

1 - 

spur r 

i 

— ipm t 

1 

— efth — r 

i 

TNLV 


163 

r t™" 
1 
1 

85,46 

I 

87.53 ! 

88.24 

80.47 


200 

1 

[ 

100.00 

100.68 

101.38 

93.03 


250 

» 

1 

95.87 

102.19 

102.89 

94.33 


315 

I 

1 

95.11 

98.03 

98.73 

89.67 


400 

1 

1 

95.28 

96.37 

97.07 

88.07 


500 

1 

92.83 

93.17 

93.87 

84.71 


630 

! 

1 

93.34 

89.72 

90.42 

80.58 


800 

1 

85.14 

83.32 

84.02 

73.92 


1000 

1 

80.39 

79.67 

80.37 

70.02 


1250 

1 

77.02 

75.25 

75.94 

65.35 


1600 

1 

71.62 

70.48 

71.17 

60.34 


2000 

1 

70.01 

68.78 

69.49 

58.03 


2500 

1 

1 

69.52 

68.73 | 

69.43 

57.50 


3150 

1 

1 

67.24 

64.23 i 

64.93 

53.25 


4000 

1 

64.88 

58.69 

59.39 

47.46 


5000 

! 

65.85 

58.36 

59.06 

46.61 


6300 

1 

66.12 

64.90 

65.60 

53.13 


8000 

1 

64.20 1 

60.44 1 

61.15 1 

48.23 

OVERALL 

VALUES 

(dB) 





SPLT 

IMTH 


DFTH 




104.05 

106.32 

107.02 





If an assumed value for a of 0.1 for all 1/3 octave bands 
[29] is used the agreement in sound pressure-level is within 
about 3 dB. However, if measured values of a are used the agree- 
ment is within about 5 dB . 


The disagreement between theory and measurement of the sound 
pressure level might be caused by several sources of error. One 
source might be caused by turning the noise source off and on. 
It was found that, in order to supply the electro-pnuematic 
driver with about 20 psia which was needed to bring the overall 

C - 2 - 
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level up to about 130 dB, the pressure to the driver fluctuated 
about 101. Another possible error lies in the intensity measure- 
ment. As discussed in the previous section, the sweeping on the 
stiffened panel was made over an imaginary surface about 5mm from 
the beams and the beam width + 5mm in between the beams. This 
might have lead to some contamination from other subareas for 
almost all frequencies. A third problem might be that the 
sources are to some extent coherent and therefore assumption (4) 
in section 4.2 is violated. 

In summary, the agreement between experimentally and 
theoretically determined sound pressure levels is within 3 to 5 
dB depending upon the value of used in the prediction scheme. 
We plan more experiments in an attempt to reveal the problem. 
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1, INTRODUCTION 

During recent years, there has been a growing interest in 
the interior noise levels of general aviation aircraft, The 
noiee is associated primarily with sources such as the propel!*- 
ere, the engines, and the engine exhaust. Studies of general 
aviation aircraft have indicated that the noise is transmitted 
into the interior through airborne and structure -borne paths. 
At the same time, there has been an increasing demand on the gen- 
oral aviation industry to reduce the cabin noise levels. Hence, 
it is important to study the noise transmission characteristics 
of such aircraft structures. 


The sound transmission property of a structure may be 
evaluated from a Knowledge of its transmission loss, which is 
proportional to the ratio of the incident sound power to the 
transmitted sound power. Hence development of both accurate 
prediction and evaluation techniques is very important in the 
study of sound transmission through fuselage walls. Here in the 
second part of this report, the fuselage structure has been 
idealized into a cylindrical shell. The dimensions of the 
cylinder were designed to simulate a half -scale model of a light 
aircraft fuselage. The transmission loss was theoretically 
estimated by using Statistical Energy Analysis (SEA). The param- 
eters involved in SEA model, such as modal densities and radia- 
tion efficiencies were obtained from the wavenumber diagrams. 
The theoretical prediction will be discussed in Section 2. 
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Section 3 discusses the experimental evaluation of Bound 
transmission lose through the cylinder. The transmitted inten- 
sity was measured using the two-microphone acoustic intensity 
technique. The results obtained from the experiments compared 
well with those from the theoretical predictions* 


Z , THEORETICAL MODEL - SOUND TRANSMISSION THROUGH A CYLINDRICAL 
SHELL 

2 . 1 Natural Frequencies of a Cylindrical Shell 

The vibration analysis of a cylindrical shell is one of the 
most widely investig''* 1 problems in engineering. Many past suc- 
cessful efforts in obtaining information on the natural frequen- 
cies and mode shapes of different cylinders with different boun- 
dary conditions have been summarized by Leissa [1]. Convention- 
ally, the equation of motion of a circular cylindrical shell i3 
derived from Love's equation. For example, for a cylinder which 
is simply-supported at the two ends, the general solution for the 
motion can be assumed to be: 

u (x,e) - A cos 2^ cos n <e-4>) , 

X Lj 

U e (x,e) - B sin sin n (B-0) , (2.1) 

U 3 (x,e) - C sin 5^1 cos n (Q-P) , 
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where U x ^ U 0 , U 3 are displacements in longitudinal, circumferen- 
tial and transverse directions respectively; m, n are mode 
numbers, and x,6 are coordinates. Theoretically, a solution for 
the natural frequencies can be found when the eigen functions in 
Equation 2.1 are substituted into the equation of motion. Unfor- 
tunately, the solutions are too complicated to be applied in real 
engineering problems [2J. When the transverse vibration is con- 
sidered exclusively, the formula to calculate the natural fre- 
quencies can be largely simplified using some practical assump- 
tions. This simplified equation is known as the Donnell- 
Mushtar i-Vlasov equation. 

2.1.1 Donnell-Mushtar i-Vlasov Equation 

Donnell 13,43 and Mushtar i [5] simplified the equation of 
motion for a thin cylindrical shell by focusing attention on the 
problem of transverse deflection, and made the following assump- 
tions [2], 

a. The contributions of in-plane deflection can be neglected in 
bending strain expressions, but not in the membrane strain 
expressions . 

b. The influence of inertia in the in-plane direction is negli- 
gible. 

c. The transverse shear terms can be neglected. 

d. The mode number n, the length L, and the radius R of the 
cylinder, satisfy the inequality such that [ 6 ]. 

(nW) » (x/l ] 2 

where X and C are the constants used in the solution 
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U 3 - Ce Xx/L coe [n(e-*)j 


The natural angular frequencies of a circular cylindrical 
shell can then be reduced to 


0) 

tnn 



Bhir 

t 4_4 

l» n 



+ D 



( 2 , 2 ) 


where p m is the mass density of material, h is the thickness, m 
and n are mode numbers in the longitudinal and circumferential 
direction respectively, and E and D are Young's modulus and bend- 
ing stiffness. The only variable t) which is the roots of the 
analogous beam equation, can be determined from the boundary con- 
ditions . 


Soedel [7] in his recent paper claimed that the assumption 
(d) in the previous paragraph is not absolutely necessary in most 
circular cylindrical shell cases. He modified the formula for 
the natural frequencies into the following form; 




(2.3) 


For the case of simply-supported boundary conditions at both 
ends, the above formula can be reduced to: 



» (2.4) 
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where m is Poisson's ratio. Szechenyi [8} derived the same equa- 
tion in a different form 



Or ] 4 M 

[(r ) 2 ♦ 


As discussed in Section 2.1.2, Equations 2.3, 2.4, and 2.5 
are formulae to calculate the natural frequencies of a simply- 
supported cylindrical shell. It was also mentioned that the 
undecided variable i) m of Equation 2.2 varies with different boun- 
dary conditions. The variable y is the roots of the analogous 
beam equation. Tests are given to see how the boundary condi- 
tions affect the natural frequencies. 

Figure 2.1 shows the differences in natural frequencies 
between a simply-supported cylinder and a clamped cylinder. It 
is observed that there exists a considerable difference for mode 


m 


(i-m 2 )r ; 


2.1.2 Boundary Conditions 



Fig. 2.1 

Natural Frequencies of a Simply-Supported Cylinder t — — ) and a Clamped Cylinder 
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numbers below five. However , for higher mode numbers, the 
differences decrease sharply. Hence the boundaries of the 
cylinder under investigation are assumed to be simply-supported 
even though they may be nearer to the clamped case. Hopefully, 
the error that results from this approximation is only limited to 
the lower modes . 

2.2 The Wavenumber Diagram of a Cylinder 

The wavenumber diagram (sometimes called the "k-space 
diagram") has been found to be a very powerful tool in the Sta- 
tistical Energy Analysis (SEA). Instead of solving for each mode 
individually, SEA investigates the average behavior of variables 
over a frequency band. Among the most important parameters in 
SEA are the modal density, mode classification (surface modes, 
edge modes, and corner modes), and loss factors. All these 
parameters may be obtained from the wave number diagram. One 
other parameter, namely, the damping loss factor of the structure 
may have to be obtained experimentally. 


2.2.1 Modal Density of Cylindrical Shell 

The wavenumber functions are defined as [9] : 


k a - (gSj [h 2 R 2 /X2 j 1-« 2 J ] 1/4 , 

and 


V 
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k c ■ (ill [hV/12(lV)] V4 , <2.7) 

where k ft and k c are the longitudinal wavenumber and circumferen~ 
tial wavenumber functions, respectively. The advantage to use 
the wavenumber functions rather than frequencies is the further 
simplification of Equation 2,5 into the following form: 


<j " 

mn 




k cl 


1/2 


’ 2rrf 

r 


f (2.8) 


where f f is the ring frequency of the cylinder; 

f « [JL 

*r 12 wR 

If the natural frequencies are normalized by the ring frequency, 
Equation 2,8 becomes ? 



Equation 2,9 gives the relationship between normalized frequency 
»> and the two directional wavenumbers k and k . With a given 
frequency v q every value of k Q has a corresponding value of k . 
Therefore, the frequencies can be plotted on a wavenumber diagram 
as shown in Figure 2.2. 

A frequency band, which consists of both upper and lower 
limits of frequencies, is plotted on the wavenumber diagram as a 
strip between two frequency contours. Each mode is represented 
by a point on the diagram. All the modes are approximately 




CIRCUMFERENTIAL WAVENUMBER 
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distributed as a regular lattice on the diagram, with each mode 
occupying an area equal to: 


yrh/L 

[l2(l-M 2 JJ 1/2 

Following thiB concept, the number of modes in a frequency band 
is then equal to the ratio of the area occupied by the strip con* 
taining the frequency band to that of a unit mode. The modal den- 
sity in a frequency band, given the unit as modes / Hz, can then 
be obtained from the number of modes and the bandwidth. Figure 
2.3 gives a plot of the modal density versus freruency for the 
cylindrical shell under investigation. 

2.2.2 Acoustic Wavenumbers 

On the wavenumber diagram, for each given frequency, v Q , the 
acoustic wavenumber k is plotted as a segment of a circle with a 
cadiua v 0 W £ c- An >' chan 3 e in elthsr the geometry or the 
material of the cylinder alters the ring frequency f and the 
critical coincidence frequency £ r » Consequently, the radius of 
the circle required to plot the acoustic wave number must also be 
adjusted. Figure 2.2 shows the plots of the structural 
wavenumbers and acoustic wavenumbers corresponding to the follow- 
ing three values of the normalized frequency v q , namely 
v q - 4.0, 3.2, 1.6, 1.0, and 0.4. 


zt LU 



• • 
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2.2,3 Acoustically Past and Slow Modes 

Acoustically fast modes of a cylinder are defined as the 
modes for which both the longitudinal wavenumber and the circum- 
ferential wavenumber are less than the corresponding acoustic 
wavenumber. These modes play a key role in the sound radiation 
from the cylinder. On the other hand, the acoustically slow 
modes are those modes for which either the longitudinal 
wavenumber k ft and/or the circumferential wavenumber k c are larger 
than the corresponding acoustic wavenumber k, Following this 
definition, the modes on the wavenumber diagram can be easily 
classified as either acoustically fast or acoustically slow as 
Bhown in Figure 2.4. 

The radiation efficiency, which describes the relationship 
between the radiated power and the structure vibrational velocity 
is much higher for acoustically fast modes than that of acousti- 
cally slow modes. Usually, the radiation efficiency is unity for 
the acoustically fast modes. However, for the acoustically slow 
modes, where the circumferential wavenumber k is less than the 
acoustic wavenumber k, while the longitudinal wavenumber k ft is 
greater than the acoustic wavenumber k, the edges aligned In the 
circumferential direction produce so-called edge modes. Because 
of the geometric constraints of the cylindrical shell, Manning 
and Maidanik observed that edge modes do not exist in the longi- 
tudinal direction [10] . For the same reasons, corner modes do 
not exist on a cylinder. Figure 2.4 shows a frequency band which 
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consists of soma acoustically fast modes and some acoustically 
slow modes. Figure 2.5 gives two extreme cases , wherein « in the 
first case all the modes in a frequency band are acoustically 
fast as in Figure 2.5(a) and the other case in which all the 
modes are acoustically slow as shown in Figure 2.5(b). 


2.2,4 Radiation Efficiency 


The radiation efficiency indicates how efficiently sound Is 
radiated from a structure to the acoustic field. It is defined 
as (11) 


a 


F 

p cS<V 2 > 
o 


( 2 . 10 ) 


where P is the acoustic power radiated from the structure, p Q c is 
the characteristic impedance of the air, S is the surface area, 
and <V 2 > is the spatially averaged mean-square velocity on the 
radiating surface. 


The radiation efficiency of the modes resonant in a fre- 
quency hand of interest varies as follows. 

1, is unity for a frequency band that includes solely acousti- 
cally fast modes. 

2* For a band which has both acoustically fast and slow modes, 
the radiation efficiency may be assumed to be given approxi- 
mately by the ratio of the number of acoustically fast modes 
to the total number of modes resonant in that band. 
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3. Foe a band which has exclusively acoustically slow modes, 
the radiation efficiency fot those modes (circumferential 
edge modes of the cylinder) is given by Szechenyi as [9): 


a 



wL 
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1/2' 

. M* 


/2 
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( 2 . 11 ) 


where f is the frequency, f is the critical coincidence fre- 
quency, In is the natural logarithm, and the other quantities are 
defined in earlier equations. The radiation efficiency of a 
cylindrical shell with the ring frequency at 2134 Hz and the 
critical coincidence frequency at 7452 Hz is plotted in the Fig- 
ure 2.6. The radiation efficiency in Figure 2.6 is given in 
decibels with a reference value of 1.0. 


2 . 3 The Sound Transmission Loss of a Cylinder 

The study here of sound transmission loss through a finite 
cylindrical shell starts with the simplest case, wherein only the 
cylindrical shell is considered and diffraction at the ends is 
ignored, The sound transmission loss is then defined as 

„ n -> T Incident Acoustic Power 1 /0 lov 

TL " ~ ^Transmitted Acoustic Power J * ( • ) 

The sound transmitted through a structure can be divided into two 
parts; the resonant transmission and the non-resonant 
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transmission. The non-resonant transmission dominates in the 
frequency range below the critical coincidence frequency, while 
the resonant transmission increases in importance at higher fre- 
quencies. The resonant transmission can be derived from the 
power balance equation as [13]. 


T W - - 101 °9 


8?r^c^ n ( w 0 ] R 


2 

rad 


w o Ms2 ( 2R r»d + R 
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(2.13) 


where n(u» o ) is the modal density in modes/Hz, R ^ is the average 
radiation resistance of all the modes resonant in the band under 
consideration is the center frequency (of the band), M is the 
mass of the cylinder per unit area, S is the radiating surface 
area, and R mec ^ is the mechanical resistance. The non-resonant 
vibration can be much more important at the frequencies where the 
resonant radiation properties are weak. 


An incident acoustic wave forces any panel structure to 
vibrate off-resonance with a flexural wavelength equal to the 
acoustic trace wave length. Consider a structural mode with 
resonance frequency f having the same wavenumber components as 
those of the acoustic wave when incident on the structure (ie. 
same trace wavenumber vector.) This mode can be excited as fol- 
lows [13). 

(i) if f > f r the vibration amplitude (hence the 
transmission) is mass controlled, this response is said 
to be sub-coincident. 

(ii) if f * f , the mode is excited at resonance and the 
response Is said to be coincident. 


- no 


(iii) if f < f the vibration and transmission are 
stiffness-controlled and the response is said to be 
super -co inc ident . 

Thus the non-resonant transmission for any incident acoustic 
wave is mass controlled if its frequency is greater than the 
natural frequency of the structural mode having the same coord i~ 
nate wavenumbers, and stiffness controlled if the frequency is 
less than this natural frequency# In a diffuse incident sound 
field all the structural modes having wavenumbers k and k_ such 
that k 2 + k 2 < k 2 will be excited in non-resonant vibration [13], 

c| C 

The non-resonant modes of a cylindrical shell can be classified 
into mass controlled modes and stiffness controlled modes below 
the ring frequency because of the curvature effect. This is dif- 
ferent from the case of a flat plate where the mass and stiffness 
controlled regions are easily defined, because all incident waves 
will be subcoincident for all frequencies up to the critical 
coincidence frequency. Although both represent non-resonant 
vibration, the mass controlled subcoincident waves are much supe- 
rior to the stiffness controlled super-coincident waves in sound 
transmission. Accordingly, below the ring frequency, the non- 
resonant transmission calculation takes only the mass controlled 
modes into account though in that region both mass and stiffness 
controlled modes exist. For those frequencies above the ring 
frequency, since all the non-resonant modes are mass controlled, 
the non-resonant transmission is the same as that of a flat 
panel. In summary, the non-resonant transmission losses are 
given by the following two equations for the two frequency 
regions [13] . 
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for v Q < 1 (below ring frequency) 


(2.15) 
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for v Q > 1, (above ring frequency) 

The total transmission loss of diffuse incident waves on the 
cylinder is the sum (in an energy sense) of the resonant and the 
non-resonant transmission loss. The transmission loss obtained 
theoretically in decibels versus frequency is shown in Figure 

3.14, 

2 .4 Computer Programs and Discussion of the Theoretical Approach 

Previous studies on the transmission of sound through 
cylindrical shells such as the work of Manning et al [10] and 
Szechenyi [8,9] utilized some approximations in obtaining the 
wavenumber diagrams. For example Szechenyi measured the areas in 
the wavenumber diagrams required for the evaluation of modal 
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densities with a planimeter and with several approximations. 
Whereas in the work reported here a series of computer programs 
have been developed to calculate the modal densities and the 
radiation efficiencies etc, from the wavenumber diagrams. These 
results are believed to be more accurate since the approximations 
made in earlier work have been avoided. 

From Equation 2.9, for each given frequency v q , the locus v Q 
is plotted on the wavenumber diagram with k and k_ as the axis 
and ordinate of the diagram respectively. A computer program has 
been developed to solve the following equation which is an 
expanded form of Equation 2.9: 

k a + ( 4k cK + ( 6k J + + l) k a < 2 ‘ 1G > 


+ 



2v k 2 
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v 2 k 4 l 
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Three hundred values of k and the same number of values of k 

c a 

are calculated for every value of v Q , In addition, the increment 

of v has been chosen to be one-sixth of an octave. This results 
o 

in a total of forty three different values of v Q in the range 
from 0.05 to 6.4. 


The surface area under a curve on the wavenumber diagram has 
been estimated by numerical integration. Whereas these areas 
have only been measured approximately in earlier work. Hence it 
has been possible to compute accurately these areas in the 
present work. Therefore, possible computation errors were minim- 


ized . 


- 113 - 


a. Szechenyi [9] made the assumption that R^/R^ch << 1 in 
the evaluation of the non-resonant transmission loss. Our 
experiments (described in the next section) suggest that 
this assumption was poor for some cylindrical shells. Hence 
the results in the present work were computed by taking into 


consideration the 

values of 

R rad 

and 

R mech 4 

b. Again in reference 

9, 

the 

ratio 

of 

the number of mass- 

controlled modes 

to 

that 

of 

the 

total number of modes 


(required in the evaluation of modal density and radiation 
efficiency of the cylinder) has been evaluated by approxi" 
mating the areas with the arcsine term in Equation 2.14. 
Whereas the ratio of the two wavenumber diagram areas have 
been evaluated by a numerical integration scheme through a 
computer program in the work reported here. 

c. The calculation of radiation efficiency for those frequency 
bands containing some acoustically fast modes was carried 
out with computer programs without any approximation, 
Szechenyi [9] approximated the area of a frequency band in 
the wave number diagram to that of a triangle as shown in 
Figure 2.7. This assumption may cause larger errors for 
those cylindrical shells with a more curved shape. 
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3. experimental investigations 

A cylindrical shell was designed to simulate a half scale 
model of a light aircraft fuselage. In Phase I of the study, the 
following series of experiments were conducted on the cylindrical 
shell assembly. 

a. Measurement of the natural frequencies and modal density of 
the cylindrical shell. 

b. Measurement of the internal loss factors of the cylinder. 

c. Measurement of the interior sound intensity generated by a 
loudspeaker mounted inside the cylinder. 

d. Measurement of the radiation efficiency of the cylinder. 

e. Measurement of the sound intensity which is transmitted from 
the cylinder. 

3 . 1 The Cylindrical Shell Assembly 

The 0.762 m (30 inch) diameter cylindrical shell was built 
by rolling 1.6 mm (1/16 inch) thick sheet steel into circular 
shape and then welding it at the joint. Its length is 1.676 m 
(66 inches). The cylindrical shell has a pair of steel rings at 
both ends. Both the inner ring and the outer ring are 9.5 mm 
(3/8 inch) in width. After the cylindrical shell and the rings 
were assembled, the outer rings were fixed on two end plates. 
The end plates were made of 6.35 mm (1/4 inch) thick steel. This 
whole assembly as shown in the Figure 3.1 was suspended from the 
ceiling in the anechoic chamber with two chains during the exper- 


iments . 


LOUDSPEAKER 



3.1 

Cylinder Assemble 




1X7 


3,2 The Natural Frequencies 

The natural frequencies of the cylinder were measured by an 
accelerometer in conjunction with a force transducer. The 
natural frequencies were obtained by observing the peaks of the 
mobility curve and the zero crossing of the phase plot of the 
curve simultaneously. The mobility is the transfer functions 
between input force and response velocity. Both the impact 
method (using the hammer) and the impedance head method (using a 
shaker) were utilized in the measurements. Table 3.1 shows a 
comparison of the first ten natural frequencies obtained theoret- 
ically and those obtained by the two experimental approaches. 
The two experimental results compared very well with each other 
while the theoretical predictions (Equations 2.4) were different 
by a few hertz. However, as far as the modal density is con- 
cerned, the theory and the measurements matched very well. Fig- 
ure 3.2 shows the mobility curve acquired from the impact method. 
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Table 3,1 

FIRST TEN NATURAL FREQUENCIES 


Mode 

Number 

Calculated 

Measured Frequency (Hz) 

m 

n 

Frequency (Hz) 

Impact 

Impedance Head 


1 


Method 

Method 

i 

4 

79.9 

75.8 

74.4 

i 

5 

81.0 

77.4 

77.6 

i 

6 

102.8 

94.8 

94.6 

i 

3 ! 

118.0 

109.0 

110.6 

i 

7 1 

135,2 

124,8 

124.6 

2 

6 1 

154.2 

151.4 

151.2 

2 

7 1 

162.1 

161.2 

161,2 

1 

8 ! 

174.6 

172.0 

169.6 

2 

5 ! 

177.7 

177.6 

177.8 

2 

8 1 

189.8 

179.4 

179.4 


3 .3 The Damping Loss Factor 


A knowledge of the damping loss factor is most essential in 
SEA theory, Three related parameters are widely used to indicate 
the damping character of structures. They are the damping loss 
factor *? mec h the reverberation time T R , and the critical damping 
ratio C (143. The relationships between these measures are: 


C « (l/ 2 )7? mec h ' 


(3.1) 


and 


T 

R 


2,2/fy 


mech’ 


(3.2) 
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3.3,1 Difficult lea In Damping Measurement a 

Damping measurement oa« ba classified into two main 
approaches, (1) the ateady atate method and (2) the decay cate 
method. The ateady atate method uaea the concept of energy flow 
balance. The ateady atate method was not used in this project. 
The decay rate method, which waa adopted in thia investigation 
employe a level recorder to record the decay in vibration of the 
structure when the exciting force is suddenly shut off, Because 
this is a widely used method, no details of the measurement pro- 
cedure will be explained here. However, the difficulties 
Involved in measuring the decay rate are listed as following: 

a. It is quite subjective as to where to draw a straight line 
on the recorded paper strip, since the recording usually 
consists of a considerable number of "ripples. * 

b. For an extremely highly damped structure, the paper speed 
must be extremely high if the decay rate is to be obtained 
accurately. 

c. The experiments have to be repeated several times in order 
to get repetitive measurements. This is a time consuming 
process and hence a major disadvantage. 


3.3.2 Damping Measurements Using the FFT 

By applying similar ideas, the conventional decay rate 
method can be very much improved by using a high sampling rate 
A/D converter. The HP5451C Fast Fourier Analyzer has sampling 
rate available as frequently as every five micro-seconds. In 
addition, data can be processed directly on the machine by using 
Fortran programs. These features make the FFT a very powerful 
tool in decay rate measurements. 


The acceleration decay curve from an accelerometer attached 
to the structure can be expressed mathematically as [15] 

A(t)-Ct+b+n(t), (3.3) 

where A(t) is the acceleration level in decibels, t is the time. 
C and b are constants, n(t) is the "noise" signal representing 
the flutter in the decay. The decay rate C, can be determined by 
estimating (C) in such a way that n(t) becomes a noise component 
with zero mean. Hence the integral of n(t)dt with respect to t 
is eliminated. Equation 3.3 is thus reduced to 

$ C 2 

/ A(t)dt - 7 tS bt . (3.4) 

o 

Both Figures 3.3 and 3.4 illustrate the procedures explained 
above . 

Based on the assumption that the acceleration level decays 
exponentially, the decay rate in decibels/sec is the constant C 
in Equation 3.3. Therefore, the decay rate can be obtained 
either by choosing any two points on the Figure 3.4 or by fitting 
the curve with least square error to Equation 3.4. On the FFT, 
searching for the best estimate of the constant C is carried out 
by a Fortran program. 
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The Result of Integrating the Acceleration with Respect to Time. 



This improved method of measurement of the decay rate has 
proved reliable and time-saving. This method has been used to 
verify the reverberation time measurements of the reverberation 
room in the Herrick Laborator ies . The results obtained by this 
new method compared well with those obtained by the conventional 
method in most frequency ranges [16]. There seems to be some 
slight disagreement in the high frequency range between these two 
types of measurements [16]. 


3.3.3 Damping Loss Factor of the Cylinder 

Figure 3.5 is the schematic diagram of the instrumentation 
used for damping measurements. The white noise signal is gen- 
erated by a B&K 1024 Sine-Random Signal Generator, and then fed 
through a B&K 2112, One-Third Octave Filter. This signal ampli- 
fied by a Phase Linear Amplifier is supplied to a shaker to drive 
the cylinder at a point. The output acceleration signal is 
picked up by a B&K 4344 accelerometer. After amplification, this 
signal is fed through a B&K One-Third Octave Filter to the A/D 
converter of the FFT. A 6V battery provides an external trigger- 
ing source to the FFT through a double switch. Whenever the 
noise source is cut off, by means of the switch the FFT starts to 
trigger. 


The damping loss factors of the cylindrical shell under 




Fig* 3.5 

Schematic of the Experimental Set Up for the Damping 
Measurement. 

f 













investigation were measured in each one-third octave frequency 
band. Table 3.2 shows the results* 


Table 3 . 2 

THE DAMPING LOSS PACTORS OF THE CYLINDRICAL SHELL 


Center Frequency 
(Hz) 

Decay Rate 
(dB/sec) 

Damping Loss 
Factor 

100 

12.71 

0.00466 

125 

12.38 

0.00363 

160 

11.98 

0.00274 

200 

12.58 

0.00230 

250 1 

9.15 

0.00134 

315 

11.25 

0.00131 

400 

12.69 

0.00116 

500 

14.59 

0,00107 

630 

14.13 

0,00082 

800 

19.69 

0.00090 

1000 

33.87 

0,00124 

1250 

38.47 

0.00113 

1600 

36.79 

0.00084 

2000 

71.15 

0.00130 

2500 

66.78 

0.00098 

3150 

52.14 

0.00061 

4000 

50.49 

0.00046 

5000 

67.73 

0.00050 

6300 

74.67 

0.00043 

8000 

184.54 

0.00084 

10000 

373.53 

0.00137 

12500 

427.11 

0.00125 

16000 

391.47 

0.00090 


3.3.4 The Radiation Losses 


The loss factor measured using the decay technique discussed 
in the previous section is the total loss factor and includes the 
damping losses as well as the radiation losses of the structure. 
Hence the true internal loss factor should be obtained by sub- 
tracting the radiation loss factor from the total loss factors. 
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In most panel studies, the radiation losses turn out to be 
very low compared to the total damping losses in the range below 
the critical frequency [17,19,20]. This is mainly because the 
panel is heavily damped at the boundaries for either clamped or 
simply supported conditions. However, an examination of the 
total loss factors, as given in Table 3.2, leads to the conclu- 
sion that the cylindrical shell under investigation is not a 
heavily damped structure. Therefore, the radiation loss should 
be subtracted from the total loss factor so as to obtain the 
correct internal loss factor. 

The radiation loss factor i? rad is defined as 


p co 

_ . ° 
^rad *W * 


(3.5) 


where o is the radiation efficiency given by Equation 2.10, and M 

2 

is the surface density in kg/m . The measurement of the radia~ 
tion efficiency is discussed in Section 3.4.2. 
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Table 3.3 

INTERNAL LOSS FACTORS OF THE CYLINDER 


Center 

Frequency (Hz) 

Total Loss 
Factor 

Radiation 
Loss Factor 

Internal Loss 
Factor 

100 

0.00466 


0.00466 

125 

0.00363 


0*00363 

160 

0.00274 


0. 00274 

200 

0.00230 


0.00230 

250 

0.00134 

0.000420 

0*00092 

315 

0.00131 

0*000840 

0.00047 

400 

0.00116 

0.000332 

0.00083 

500 

0.00107 

0.000265 

0.00080 

630 

0.00082 

0.000337 

• 0.00048 

800 

0.00090 

0.000265 

0.00063 

1000 

0.00124 

0.000335 

0.00090 

1250 

0.00113 

0.000238 

0.00089 

1600 

0.00084 

0.000166 

0.00067 

2000 

0.00130 

0,000419 

0.00088 

2500 

0.00098 

0.000134 

0.00085 

3150 

0.00061 

0.000030 

0.00031 

4000 

0.00046 

0.0000167 

0.00044 

5000 

0.00050 

0.0000127 

0.00049 

6300 

0.00043 

0.0000150 

0.00042 

8000 

0.00084 

0.0000105 

0.00082 

10000 

0.00137 

0.0000042 

0.00013 


3.4 Sound Intensity Radiated From the Cylinder 


A 15 watt, double coned loudspeaker was mounted, on brackets 
six inches away from the endplate, at one end of the cylinder. 
It was oriented toward the end plate and was mounted eccentric to 
it. After a steady sound field was produced inside the cylinder, 
the sound intensity transmitted through the cylindrical shell was 
measured by sweeping the combination of two microphones over the 
surface of the cylinder. The intensity measurement is shown in 
Figure 3.6 (a). The instrumentation used for the transmission 
loss measurements is shown in Figure 3.6 (b) * 






Fiq. 3.6b Instrumentations for the Acoustic 

Intensity Measurements 
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3.4.1 Two Microphone Intensity Measurement 

The use of the combination of two microphones for the meas- 
urement of acoustic intensity has been discussed in Part 1 of 
this report. In this study, three different microphone spacings, 
6.35 mm, 12.7 mm, and 25,4 mm, were employed in order to produce 
satisfactory results at both high and low frequencies. For all 
the measurements, the signal-to-noise ratios were carefully 
checked. These ratios varied from 25 dfi at low frequencies to 
more than 50 dB at high frequencies. 


3,4,2 Radiation Efficiency 

As discussed in the section 2.2.4, the radiation efficiency, 
can be determined by measuring the velocity of the vibrating 
structure and the power radiated by it. In this investigation, 
the spatially averaged velocity is taken at. 60 randomly distri- 
buted locations over the cylindrical shell. Simultaneously the 
sound power radiated was measured with the two-microphone inten- 
sity method. The results of the radiation efficiency were 
analyzed in both narrow and one-third octave bands. Figure 3.7 
presents the narrow band analysis over a frequency range of 0 Hz 
to 10000 Hz. Figure 3.8 exhibits a comparison between the 
theoretical predictions (Sec. 2.2.4) and the experimental results 
of the radiation efficiencies for one-third octave center fre- 
quencies. The agreement seems to be fairly good. 
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3 * 5 Sound Transmission Through the Cylindrical Shell 

As explained in the project proposal [18), attention ie 
focused in this project on techniques to determine the sound 
transmission through an aircraft fuselage. In phase 1 of the 
study, emphasis is placed on the sound transmission through a 
cylindrical shell. 

Traditionally, the transmission loss of a structure is determined 
from the noise reduction with associated corrections for the area 
of the structure and the absorption in the receiving space. In 
order to evaluate the transmission loss, usually some complicated 
testing procedures and mathematical calculations are needed. 
Moreover, the basic assumption of two reverberant fields in the 
measurement, required either a pair of costly reverberation rooms 
or sacrifice of measurement accuracy. Now, the measurement of 
transmission loss is completely simplified by the two- microphone 
intensity method [12). 

The sound field Inside the cylinder is generated by a 
loudspeaker placed inside the cylinder. As generally known, the 
sound pressure is distributed unevenly in the low frequency 
range. In other words, there exists some standing waves. Among 
them, the first two longitudinal modes at 100 Hz and 200 Hz and 
the first two cross modes at 265 Hz and 439 Hz tend to have high 
amplitudes at their resonance frequencies. Around these modal 
frequencies (mainly below 500 Hz), the modal densities are too 
low to build up a diffuse field. Hence the assumption that it is 
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a diffuse sound field inside the cylinder may cause a large 
discrepancy in the results in the low frequency region. In order 
to obtain a good spatially averaged sound pressure level, a bent 
microphone probe is moved both longitudinally and clrcumferen- 
tially inside the cylinder during the sound level measurements. 
Table 3.4 shows the spatially averaged sound level (from 300 sam- 
ples) and the sound level at two arbitrarily picked points inside 
the cylinder. The agreement between the three measurements at 
high frequencies suggests that the sound field is indeed quite 
diffuse. 


I 


Table 3.4 

SOUND PRESSURE LEVELS OF TWO FIXED POINTS INSIDE THE CYLINDER 
COMPARED TO THE SPATIALLY AVERAGED SOUND PRESSURE LEVEL 


Frequency (Hz) 


6 
8 

10 
12 
16 
200 
25 
31 
40 
50 
63 
800 
1000 
1250 
1600 
2000 
2500 
3150 
4000 
5000 
6300 
8000 
0000 



L p at Point A 


61.3 

67.7 

87.3 

74.9 

65.3 

78.7 

79.8 
92.5 

98.0 

94.9 

96.0 

100.3 

102.1 

104.7 

100.7 

99.9 
97.8 

103.4 

100.1 

97.1 

96.1 

95.7 

88.7 


Lp at Point B 


63.0 
75.6 

83.5 

75.9 

67.5 

87.9 

77.3 

92.8 

95.9 
97.2 

98.9 

101.7 

103.5 

106.6 

102.8 

99.9 

95.9 

103.6 

101.8 

94.9 
97.2 

96.4 

90.0 


9 

7 
6 

8 
6 
2 

5 
8 
9 
3 

6 

100.1 

104 . 4 

105.5 
102.4 

99.9 

96.9 
103.9 
101.3 

95.6 

96.9 

96.8 

90.9 


Sound Intensity Radiated to the Exterior of the Cylinder 


The sound intensity transmitted through the cylinder was 
measured by the two microphone technique on the exterior surface 
of the cylinder. It is of interest to compare the sound radiated 
from different sections of the cylinder surface area . Figure 3.9 
shows the effect of variation of loudspeakers location. Areas #1 
and #12 shown in the figure are located at two opposite ends of 


the cylinder. The loudspeaker was installed near the area #12. 





FREQUENCY (Hz) 
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Differences are found in the low frequency range whereas little 
discrepancy is noticeable at high frequency. Secondly, the edge 
effect is presented in Figure 3.10. Area #1 has one edge clamped 
by two rings while area #2 is not clamped along any edge. From 
the figure, the transmitted intensity seems to be little affected 
by boundaries. Thirdly a comparison of sections with and without 
a welding joint was made. The area #8 which consists of a weld- 
ing joint seems to radiate less power than the area #2 which is 
without the welding joint. This is illustrated in Figure 3*11* 

3.5.2 Sound Transmission Loss 

The sound transmission loss is calculated by ten times the 
logarithm of the ratio of the incident intensity to the transmit- 
ted intensity. Under an assumption that the interior sound field 

is diffuse, the incident intensity can be calculated from 
2 

p /4p Q c. The transmission loss is then obtained from these 
intensity levels. Figures 3.12 and Figure 3.13 show the 
transmission loss of the cylindrical shell in a narrow frequency 
band and in one-third octave averages respectively. 

The experimental results are compared with the theoretical 
transmission loss predictions in Figure 3.14. The theory agrees 



Fig. 3.10 
Sound Inti 
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well with the experiment for frequencies above 500 Hz, It is 
especially noticed that the two dips at the ring frequency (2134 
Hz) and at the critical coincidence frequency (7452 Hz) are 
exactly predicted. However, discrepancies between the theory and 
the measurement are more pronounced at low frequency. This might 
be explained in the following ways: 

a. As discussed in Section 3.5.1, the interior sound field is 
not fully diffuse at low frequencies. Hence the use of 
p /4p Q c to obtain the incident sound intensity is an invalid 
assumption. 

b. The interior volume of the cylinder does not have enough 
modes in a frequency band of interest for coupling with the 
structural mode in the low frequency range. The SEA theory 
can be completely wrong if none or very few resonant modes 
exist in certain frequency bands. 

c. When the spacing between two microphones is not large 
enough, the intensity method might result in intensity 
errors at very low frequencies. This is because of the fact 
that in measurements with two microphones insufficiently 
spaced, the phase information of a low frequency wave (large 
wave length) can easily be in error. 
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